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Abstract 



Wc propose a probabilistic cellular automata model for the spread of inno- 
vations, rumors, news, etc. in a social system. The local rule used in the 
model is outertotalistic, and the range of interaction can vary. When the 
range R of the rule increases, the takeover time for innovation increases and 
converges toward its mean-field value, which is almost inversely proportional 
to R when R is large. Exact solutions for i? = 1 and R — oo (mean-field) are 
presented, as well as simulation results for other values of R. The average 
local density is found to converge to a certain stationary value, which allows 
us to obtain a semi-phenomenological solution valid in the vicinity of the 
fixed point n = 1 (for large t). 



1 Introduction 



Diffusion phenomena in social systems such as spread of news, rumors or in- 
novations have been extensively studied for the past three decades by social 
scientists, geographers, economists, as well as management and marketing 
scholars. Traditionally, ordinary differential equations have been used to 
model these phenomena, beginning with the Bass model and ending with 
sophisticated models which take into account learning, risk aversion, nature 
of innovation etc. ||, |T0|, |1T|. Models incorporating space and spatial dis- 



tribution of individuals have been also proposed, although most research in 
this field has been directed to the refining of discrete Hagerstrand models 
0, P and to constructing partial differential equations similar to diffusion 
equations known to physicists [0. 

Diffusion of innovations (we will use this term in a general sense, meaning 
not only innovations, but also news, rumors, new products etc.) is usually 
defined as a process by which the innovation "is communicated through cer- 



tain channels over time among the members of a social system" |jTT|. In most 
cases, these communication channels have rather short range, i.e. in our de- 
cisions we are heavily influenced by our friends, family, coworkers, but not 
that much by unknown people in distant cities. This local nature of social 
interactions makes cellular automata a convenient tool in modeling of diffu- 
sion phenomena. In fact, epidemic models formulated in terms of automata 
networks have been succesfully constructed in recent years 0, |5, . 

In what follows, we will propose a model for the spread of innovations 
formulated in terms of probabilistic cellular automata, and investigate the 
role of the range of interaction in the diffusion process. 

Let us define a symbol set Q = {0, 1, ...k — 1} and a bisequence over Q 
which is a function on Z (the set of all integers) to Q, that is, s : Z — ^. The 
set of all bisequences over Q, i.e. S = Q^, is called the configuration space, 
and its elements are called configurations. The mapping T : S —>■ S defined 
by [T(s)](z) = s{i + 1) for every z G Z is called the (left) shift. 

Let F : S S he a continuous map commuting with T. In general, a 
cellular automaton ( CA ) [|12| is a discrete dynamical system defined by 

St+i = F{st). (1) 

As proved in , in this case one can always find a local mapping such that the 
value of St+i{i) depends only on a finite number of neighboring St{j) values. 
In the remaining part of this paper we will only consider 2-state probabilistic 
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cellular automata, with a dynamics such that St+i{i) depends on St{i) and 
(7t(z), where 



^t{i) = Yl st(i + n)p(n), 

n=—OD 

and p is a nonnegative function satisfying 

oo 



(2) 



(3) 



In the simplest version of our model, the sites of an infinite lattice are all 
occupied by individuals. The individuals are of two different types: adopters 
{st{i) = 1) and neutrals {st{i) = 0). Once an individual becomes an adopter, 
he remains an adopter, i.e. his state cannot change. At every time step, each 
neutral individual can become an adopter with a probability depending on 
the parameter a. In this paper we will assume that this probability is equal 
to a. 

The model can be viewed as a probabilistic cellular automaton with prob- 
ability distribution 



The transition probability Pb^a is defined as 

Pb^a = P{st+i{i) = b\st{i) = a), 



(4) 
(5) 



(6) 



and represents the probability for a given site of changing its state from a 
to b in one time step. The transition probability matrix in our case has the 
form 





(Ttii) 1 



(7) 



Subsequently, we will consider an uniform outertotalistic neighborhood of 
radius R, with the parameter a defined as 



-1 R 

= I Z ^t{i + n) + J2 ^t{i + n)] . 

\n=-R n=l 



(8) 



a is then a local density of adopters over 2R closest neighbors. This choice 
of cr, although somewhat simplistic, captures some essential features of a 
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real social system: the number of influential neighbors is flnite and these 
neighbors are all located within a certain flnite radius R. Opinions of all 
neighbors have equal weight here, which is maybe not realistic, but good 
enough as a flrst approximation. Let rit be the global density of adopters 
at time t (i.e. number of adopters per lattice site), and = 1 — rit. Since 
Pi^i = 1, the number of adopters increases with time, and lim^^oo = 1. 
If we start with a small initial density of randomly distributed adopters no, 
nt follows a characteristic S-shaped curve, typical to many growth processes. 
The curve becomes steeper when R increases, and if R is large enough it 
takes only a few time steps to reach a high n (e.g. n = 0.99). Figure |1| shows 
some examples of curves obtained in a computer experiment with a lattice 
size equal to 10^ and = 0.02. 



2 Exact Results 

Average densities nt+i and nit+i can be obtained from previous densities nt 
and rrit using the transition probability matrix ( () denotes a spatial average) 



nt+i 




hence: 



nt+i = nt{Pi^i) + (1 - nt){Pi^o), 
and using (|^ we have 

nt+i = nt + (1 - nt){at{i)). 



(9) 
(10) 

(11) 



This difference equation can be solved in two special cases, R = 1 and R = oo. 

If i? = 1, only three possible values of local density a are allowed: 0, | 
and 1. The initial conflguration can be viewed as clusters of ones separated 
by clusters of zeros. Only neutral sites adjacent to a cluster of ones can 
change their state, and they will become adopters with probability |. All 
other neutral sites have local density equal to zero, therefore they will remain 
neutral, as shown in an example below (sites that can change are underlined): 
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Figure 1: Density of adopters (vertical axis) as a function of time (horizontal 
axis) for several values of radius R 
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This implies that the length of a cluster of zeros will on average decrease by 
one every time step, i.e. 

M{s,t + 1) = M{s + l,t), (12) 

where M(^s,t) denotes a number of clusters of zeros of size s at time t, and 
furthermore 

M{s,t) = M{s + t,0). (13) 

For a random configuration with initial density no, the cluster density is given 
by 

M{s,0) = {l-noynl, (14) 

hence 

M{s,t) = {l-noy+'nl. (15) 
The density of zeros can be now easily computed as 

oo 

mt = J2sM{s,t) = {l-noy-^\ (16) 

s=l 

and finally 

n, = l-{l-noY^\ (17) 
Using the above expression, we obtain 

nt+i = nt + {1- nt)no, (18) 

Comparing with (|lT]) this yields {at) = uq, i.e. the average probability that 
a neutral individual adopts the innovation is time independent and equal to 
the initial density of adopters. 

When i? — > oo, local density of ones becomes equivalent to the global 
density, thus in (0) we can replace {at) by rit'. 

nt+i = nt + {1- nt)nt, (19) 

or 

1 - ?7,+i = (1 - nty. (20) 

Hence 

nt = l-il-nof. (21) 

Note that this case corresponds to the mean-field approximation, as we ne- 
glect all spatial correlations and replace the local density by the global one. 
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3 Simulations 



Solutions obtained for the limiting cases discussed in the previous section 
suggest that in general the density of adopters might have the form 

n, = l-{l-n,y^''^\ (22) 

where / is a certain unknown function satisfying /(t, 1) = t+1 and /(t, oo) = 
2*. Plots of this function for several different values of R are shown in Fig- 
ure ^. They were constructed by measuring Ut and using 

Two apparent features of /(t, R) can be observed: for every finite i?, the func- 
tion /(t, R) becomes asymptotically linear, and the slope of the asymptote 
increases with R. This slope can be defined as a limit 

a{R) = hm (24) 

and as Figure ^ indicates, a(i?) increases with R. 

Since /(t, oo) = 2*, we can "rescale" f{t,R) and introduce g{t,R) = 
2~^f(t,R), so that limn^oo g{t, R) = 1- Several sample graphs of g{t,R) as 
a function of R for different values of t are shown in Figure This set of 
graphs reveals certain regularity of g{t,R), namely the shape of the g(t,R) 
for a given t is the same as the shape of g{t + 1, R) with vertical axis rescaled 
by a factor 2, i.e. g{t + 1, 2R) = g{t, R). In terms of our original function / 
this becomes 

/(t + fc,2'=/?) = 2'=/(t,/?), (25) 

where k is a positive integer. This condition is more accurate for larger R, 
although we found it well satisfied even for rather small R (less than 10). It 
is seriously violated only for R = 1, when f(t,l) = t + l. 

Combining ([TT]) and (P^, we can express average local density a in terms 
of fit, R): 

(^^) = l_(l_^,)/(m,ii)-/(<,ii)_ (26) 

Since f{t, R) is linear for large t, the average local density must become 
time-independent for large t, and therefore 

lim(ai) = !-(!- no)"(^). (27) 
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Figure 2: Plots of f{t, R) (vertical axis) as a function of time (horizontal 
axis) for several values of radius R. 
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Figure 3: Plots of g{t, R) = 2 ^f{t, R) (vertical axis) as a function of radius 
R (horizontal axis) for several values of time t. 
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Figure 4: Typical plot of the at versus time, for R = 15. 

This limit will be referred to as (o"oo)- The average local density converges 
to ((Too) rather quickly. For example, when R = 15, only 10 time steps are 
needed in practice to reach the limit value, as shown in Figure ^. Condition 
(P^l ) can be now used to find the dependence of ctqc upon the radius R. As 
we already noted, f{t, R) is asymptotically linear when t —>■ oo, and the 
convergence is so fast that we can assume f{t, R) = a{R)t + b{R) for large t. 
Using condition (^) we have 

f{t, 2R) = 2f{t -l,R) = 2a{R)t + b{R) - a{R), (28) 

which implies that a{2R) = 2a{R), i.e. a{R) is proportional to R. This can 
be easily verified if we rewrite (|27| ) as 

a(i?)ln(l-no)=ln(l-(aoo)), (29) 

and plot the right hand side of the above equation versus R, as shown in 
Figure ^ Of course, the slope a{R) could also be measured directly from 
graphs of f(t, R), and then plotted as a function of R. In both cases, a{R) = 
jR, where 7 = 0.89 ± 0.05, and finally 

(aoo) = 1 - (1 - no)^^. (30) 
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Figure 5: Plot of ln(l 



((Too)) as a function of R. 



As before, this expression is not very accurate if the radius R is small (close 
to 1). 

The value of a{R), in a sense, provides us with a measure of a "speed" of 
convergence toward fixed point n = 1, assuming we are sufficiently close to 
n — 1. It does not carry, however, useful information about the early stage 
of the process, away from n — 1 (all our previous derivations assume large 



Therefore, we will consider another quantity, often used in marketing, 
called the takeover time, and usually defined as the time required to go from 
n = 0.1 to n = 0.9. For an S-shaped curve, the inverse of its maximal slope 
carries a similar information as the takeover time, therefore we define 



Subsequently, Tr will be used as a measure of the takeover time for a model 
with neighborhood radius R, and for simplicity we will just call it takeover 
time. 



t). 




(31) 
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Figure 6: Asymptotic behavior of Tr — T^o for large R. 

For R = oo, the largest n^+i — Ut occurs at t = 5, therefore 

Too = K - n,)-' = ((1 - nof' - (1 - nof'Y' (32) 

We used uq = 0.02, which yields Too = 4.00915. 

We studied the asymptotic behavior of Tr — T^o as R tends to infinity 
(Figure H). This log- log plot shows that Tr tends to Too as where 
a = 1.03 ± 0.01. This exponent depends rather weakly on the initial density 
of adopters no. For example, if uq is two times smaller (no = 0.01), a changes 
its value by approximately 2%, i.e. a = 1.05 ± 0.01. Note that in this case 
Too is different too, namely Too = 4.01051. 

4 Conclusion and Remarks 

We have studied a probabilistic cellular automata model for the spread of 
innovations. Our results emphasize the importance of the range of interaction 
between individuals, as the innovation propagates faster when the radius is 
larger. Increased connectivity between individuals reduces constraints on 
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the exchange of information, and thus the growth rate increases. It should 
be pointed out that larger connectivity can be also achieved by increasing 
dimensionality of space. We performed some measurements of the takeover 
time in two and three dimensions (to be reported elsewhere), and as expected, 
the takeover time was found smaller when dimensionality increased. This is 
also consistent with the fact that, in general, in higher dimensions we are 
closer and closer to the mean-field approximation. 

The behavior of local and global densities of adopters is an interesting 
feature of the model. While the global density of adopters always converges 
to 1, the local density does not. It quickly reaches certain "equilibrium value" 
and stays constant after that (for infinite lattice only, since for a finite lattice 
both densities reach 1 after a finite number of steps). 

More sophisticated (and realistic) models can be constructed by relaxing 
some assumptions discussed in the introduction. For example, since every 
technology has a finite life span, we can allow adopters to go back to neutral 
state with a certain probabihty. Moreover, we can incorporate "reluctance" 
to adopt by assuming that the probability of adoption is not equal, but pro- 
portional to the local density of adopters. Such a generalized model exhibits 
a second order phase transition with a transition point strongly dependent 
on R. These results will be reported in details elsewhere. 
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